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We formulate a Lippmann-Schwinger- type resonating-group equation to calculate invari- 
ant amplitudes of the quark-model baryon-baryon interaction. When applied to our recent 
SUe quark model for the nucleon-nucleon and hyperon-nucleon interactions, this technique 
yields very accurate phase-shift parameters for all partial waves up to the energies of several 
GeV. The technique also has a merit of a straightforward extension to the G-matrix equa- 
tion. A new analytic method is proposed to calculate the quark-exchange Born kernel for 
the momentum-dependent two-body interaction. The partial-wave decomposition in the mo- 
mentum representation is carried out numerically. The invariant amplitudes are then used 
to calculate single-nucleon potentials in normal nuclear matter for high incident momenta 
qi > 3 fin -1 , in which the so-called t e ^ p prescription is found to be a good approximation 
to the single-particle potentials directly calculated in the lowest-order Brueckner theory. 

§1. Introduction 



O ' Though the quantum chrome-dynamics (QCD) is believed to be the fundamental 

theory of the strong interaction, it is still too difficult to apply it directly to two- 
baryon systems. At this stage a number of effective models have been proposed to 
understand the nucleon-nucleon (NN) and hyperon-nucleon (YN) interactions from 
^ ■ basic elements of quarks and gluons.EP Among them the non-relativistic quark model 

has a unique feature that it enables us to take full account of a dynamical motion 
of the two composite baryons within a framework of the resonating-group method 
(RGM).EP The model describes confinement with a phenomenological potential and 
uses quark-quark (qq) residual interactions consisting of a color analogue of the 
Fermi-Breit (FB) interaction. In the last several years, it was found that a proper 
incorporation of the meson-exchange effect is essential to make such a model realistic 
for the description of the NN and YN interactions. &.0)»&»B»0).i>.&.@) 

We have recently achieved a simultaneous description of the AW and YN in- 
teractions in the RGM formulation of the spin-flavor SUq quark model. lP' 
In this model the meson-exchange effect of scalar (S) and pseudo-scalar (PS) meson 
nonets is incorporated in the quark Hamiltonian in the form of effective meson- 
exchange potentials (EMEP) acting between quarks. The flavor symmetry breaking 
for the YN system is explicitly introduced through the quark-mass dependence of 
the Hamiltonian, as well as the flavor dependence of the exchanged meson masses. 
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An advantage of introducing the EMEP at the quark level lies in the stringent re- 
lationship of the flavor dependence, which is revealed in the various pieces of the 
NN and YN interactions. In this way we can utilize our rich knowledge of the NN 
interaction to minimize the ambiguity of model parameters, which originates from 
the scarcity of the present experimental data for the YN interaction. 

We already have three different versions called RGM-FeP'O 1 , FSSB'^-tP and 
RGM-H §' El, differing in the treatment of the EMEP. The model called RGM-F 
introduces, besides the central force of the S-meson nonet, only the tensor component 
generated from the n- and iT-meson exchanges, and uses some approximations in 
evaluating the spin-flavor factors of the quark-exchange RGM kernel. On the other 
hand, FSS and RGM-H calculate the spin-flavor factors explicitly at the quark level, 
and include the spin-spin terms originating from all members of the PS-meson nonet. 
The SU3 relation of the coupling constants emerges as a natural consequence of the 
SUq quark model. For S-mesons, the F j (F+D) ratio turns out to take the SUq value 
of purely electric type. This is too restrictive to reproduce existing experimental data 
for the low-energy YN cross sections. We relax this restriction in two ways; one is 
to change the mixing angle of the flavor-singlet and octet scalar mesons only for 
the UN (I = 3/2) channel, and the other is to employ the same approximation as 
RGM-F solely for the isoscalar S-mesons, e and S* . We call these models FSS and 
RGM-H, respectively. Predictions of these two models are not very different except 
for the roles of the LS^ force in the AN - EN (I = 1/2) coupled-channel system. 
The SU3 parameters of EMEP, S-meson masses, and the quark-model parameters are 
determined to fit the NN S-wave and P-wave phase shifts under the constraint that 
the deuteron binding energy and the 1 5o scattering length are properly reproduced. 
The low-energy cross section data for YN scattering are also employed to fix the 
parameters in the strangeness sector, especially the up-down to strange quark-mass 
ratio A = m s /m u d- The reader is referred to the original paperstP'Q' for a full account 
of the models and the model parameters. 

So far we have solved the coupled-channel (CC) RGM equation in the improved 
variational method developed by KamimuraEJ). In this method each Gaussian ba- 
sis function is smoothly connected to the positive- or negative-energy asymptotic 
waves, which are obtained by numerically solving a "local" CC Schrodinger equation 
consisting of the long-range one-pion tensor force and the other EMEP. Although 
this technique gives accurate results at laboratory energies up to about 300 MeV, it 
seems almost inaccessible to higher energies due to the rapid oscillation of the rela- 
tive wave functions. In this paper we formulate an alternative method to solve the 
CC RGM equation in the momentum representation; namely, we derive a Lippmann- 
Schwinger-type RGM equation which we call an LS-RGM equationH In this method 
all the necessary Born amplitudes (or the Born kernel) for the quark-exchange kernel 
are analytically derived by using a new transformation formula, which is specifically 
developed for momentum-dependent two-body interactions acting between quarks. 
The partial wave decomposition of the Born kernel is carried out numerically in the 



*' The idea to solve the RGM equation in the momentum representation, in order to avoid the 
rapid oscillation of relative wave functions at higher energies, is not new. See, for example, Ref.[L2). 
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Gauss-Legendre integral quadrature^JThe LS-RGM eauation is then solved by using 
the techniques developed by NoyesE^P and KowalskitS. Although this method re- 
quires more CPU time than the variational method, it gives very stable and accurate 
results over a wide energy range. Since we first calculate the Born amplitudes of the 
RGM kernel, it is almost straightforward to proceed to the G-matrix calculation. 

As an application of the present formalism, we discuss single-particle (s.p.) po- 
tentials in normal nuclear matter. This application is motivated by the G-matrix 
calculation of the AW and YN interactions in the continuous prescription of inter- 
mediate spectral!^. The s.p. potentials of the nucleon and hyperons predicted by 
the model FSS have a flaw that they are too attractive in the momentum region 
qi = 5 ~ 20 fm^ 1 . We will show that this particular feature of FSS is related to 
the ill-behavior of the spin-independent central invariant amplitude at the forward 
angles, which is traced back to too simple S-meson exchange EMEP in this model. 
We analyze this problem by using an approximation of the s.p. potential in the 
asymptotic momentum region, in terms of the T-matrix solution of the LS-RGM 
equation. This technique is sometimes called the t e ^p prescription. We find that 
this procedure gives a fairly good approximation of the s.p. potentials predicted 
by the G-matrix calculation, even for such a small momentum as q± ~ 3 fm -1 
(Tiab ~ 200 MeV), as long as the real part of the s.p. potential is concerned. 

In the next section we start from the standard RGM equation for the (3q)-(3q) 
system and derive the LS-RGM equation in the momentum space. Appendix A gives 
a convenient formula to calculate the Born kernel of the quark-exchange RGM kernel. 
The formula is especially useful in the case of the momentum-dependent qq force. 
The Born kernel of FSS is explicitly given in § 2.2 and in Appendix B. Since the 
flavor dependence of the spin-flavor factors plays an essential role in the application 
of the present formalism to the AW and YN interactions, the structure of the Born 
kernel and the scattering amplitudes is carefully described in § 2.3 in terms of the 
Pauli-spinor invariants. Some useful expressions of the partial wave decomposition 
of the Born kernel and the invariant amplitudes are given in Appendices C and D, 
respectively. The t e ^ p prescription is derived in § 2.4 as a method to calculate the 
asymptotic behavior of the s.p. potentials in the G-matrix calculation. In § 3 we 
compare the phase-shift parameters obtained by the LS-RGM formalism with those 
by the improved variational method. The system we choose as an example is the 
most complicated system of the AN - UN (I = 1/2) channel coupling. The invariant 
amplitudes and the scattering observables of the AW system in the intermediate 
energies between T/ a & = 400 ~ 800 MeV are discussed in § 4. The s.p. potential of the 
nucleon in normal nuclear matter is calculated in § 5, by using the t e ^ p prescription. 
It is shown that the behavior of the spin-independent central invariant amplitude at 
the forward angles is related to the attractive behavior of the s.p. potentials around 
qi = 5 ~ 20 fm -1 region. Also shown is a preliminary result of s.p. potentials in 
an improved model, which incorporates momentum-dependent higher-order terms of 
the S and vector (V) meson EMEP. Finally, § 6 is devoted to a brief summary. 
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§2. Formulation 

2.1. LS-RGM equation 

As mentioned in the Introduction, FSsll'&'l^, RGM-H&'H and their preceding 
version RGM-FEreP are formulated in the (3q)-(3q) RGM applied to the system of 
two (Os) 3 clusters. The qq interaction is composed of the full FB interaction with 
explicit quark-mass dependence, a simple confinement potential of quadratic power 
law and EMEP acting between quarks. The RGM equation for the parity-projected 
relative wave function Xa(R-) is derived from the variational principle {S^E— H\&) = 
0, and it reads as& 



£a + 



h 2 



(-)' 

2fM a \dRj 



Xl(R) =T,J dR ' G aa ,(R,R';E) xl'(R') , (2-1) 



where G aa >(R, R'; E) is composed of various pieces of the interaction kernels as well 
as the direct potentials of EMEP: 



G aa / {R, R j E) — 5(R — R ) 



3 & 



Q 



M ( a n J(R,R!)-e Q M% a ,(R,R') . (2-2) 



The subscript a stands for a set of quantum numbers of the channel wave function; 
a = [l/2(ll)oi,l/2(ll)o2] SS Z YII Z ;V, where l/2(ll)a is the spin and SU 3 quan- 
tum number in the Elliott notation (A//), a{= YI) is the flavor label of the octet 
baryons (JV = Ul/2), A = 00, U = 01 and 3 = -1(1/2)), and V is the flavor- 
exchange phase. W In the NN system with a\a2 = NN, V becomes redundant since 
it is uniquely determined by the isospin as V = (-l) w . The relative energy e a m 
the channel a is related to the total energy E of the system through e a = E — E™ 1 
with E™ 1 = E™* + In Eq. (2-2) the summation over Q for the exchange ker- 

nel -M- f^i involves not only the exchange kinetic-energy (K) term, but also various 
pieces of the FB interaction, as well as several components of EMEP. The FB in- 
teraction involves the color-Coulombic (CC) piece, the momentum-dependent Breit 
retardation (MC) piece, the color-magnetic (GC) piece, the symmetric LS (sLS) 
piece, the antisymmetric LS (aLS) piece, and the tensor (T) piece, The EMEP 
contains the central (CN) component from the S-mesons, and the spin-spin (SN) 
and tensor (TN) terms originating from the PS mesons. The contribution from a 
particular meson exchange is denoted by (3 in Eq. (2-2) only for the direct part of 
G aa >(R,R';E). The explicit form of these qq forces is given in Refs.|6|) and 16), and 
the corresponding basic Born kernel defined through 

M* a ,(q f , q t ; E) = ( e^ R \ G aa ,{R, R';E) \ e^ R ') 

= ( e iqrR Va F I G(R, R';E)\ e^'rj^F) (2-3) 
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is given in the next subsection and Appendix B. Here rj^ F is the spin-flavor wave 
function at the baryon level, which is defined in Eq. (2.9) of Ref. 17). 

To start with, we use the well-known Green function to convert the RGM equa- 
tion Eq. fl2-ip to an integral equation which has a parity-projected incident plane 
wave in the a channel: 



X-ya C^j ^(1, 



s. 



7, a 



+ {-l) Sa V a e 



1 2/x 7 



J dR ' \R-R'\ ^ J dR " ^(^-R -> E ) X*p a {R",k a ) . (2-4) 



Here it = (—l) Sa V a = (—l) s i 3 V ! 3 = (— l) s ~tV~,, because of the parity conservation. 
The asymptotic wave of Eq. (2-4) with R — > oo is given by 



■)(Z a {R, k a ) ~ 5, 



7,a 



+ (-lf a V a e 



+ - 



1 



Vj R 



M ia (k 1: k a ;E) + (-lf a V a M 7Q (fc 7 , —k a ;E) 



(2-5) 



with fc 7 = kjR and 



4^2 



Y,j dR j dRl e ~ ikl ' R G -yp( R ' R '> £ ) fc «) • ( 2 -6) 



The Born amplitude is obtained by approximating 



(2-7) 



in Eq. (2-6). Though k Q and k 7 are related to the total energy E by the on-shell 
condition 

f> 2 f? 



E = ET + 



-ki 



"> 1 2^ 



ET + 



-ki 



2-? 



it is convenient to relax this condition in order to define a more general Born am- 
plitude. Namely, fe 7 and k a are denoted by k and fc', respectively, in what follows. 
Then the Born amplitude reads asta 



M^(k,k';E) 



yj fj/syfj/Q k k 

4-Kfi 2 



M* a (k, k'; E) + (-l) s «V a M* a (k, -k'; E) 



(2-9) 



The Born amplitude Eq. ( |2 9D has a high degree of symmetries which can be directly 
derived from those of Eq. (2-3). First M® a orn (k,k';E) is a function of k 2 , k' 2 and 
cos# = k ■ k' with real coefficients, and satisfies the symmetry 

M Born (fc, k'; E) = M Born (fc, k'; Ef = M Bom (fc', k; E) . (2-10) 

*-* In Eqs. ( ]2-9| ), (2-19) and (2-21), the relative energy e 7 in the basic Born kernel (in the prior 
form) should be calculated from the total energy E through e 7 = E — Ey . 
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The parity conservation implies 



M, 



Born 



(k,k';E)=M^ a oin (-k,-k';E) 



■ya 



= (-l) 5 ^ 7 M^ a OTn (-k,k';E) 
= (-l) s -V a M^ a oin (k,-k';E) , 

with (-1) 5 ^P 7 = {-l) Sa V a = 7T. 

We now move to the momentum representation in Eq. (2-4) through 

Xjai-R) ^a) = / d-k e ^ X-yafei k a ) . 
Then Eq. (2-6) is expressed as 



M-y a (kj, k a ; E) 



(2-11) 



(2-12) 



(2-13) 



where the Born kernel V^p(k,k';E) defined by 

Vtfik, k'; E) = \[ M^(k, k'; E) + {-l)^Vp M^(k, -k'; E) 



is related to the Born amplitude through 



M^(k,k';E) 



2^ 



V 7a (k,k';E) 



We write Eq. ( gig) as 

M^ Q (k^,k a ;E) = 
and define the on-shell T-matrix by 



2irn 



(2-14) 



(2-15) 



(2-16) 



T 7a (fc 7 , k a ;E) = ^J2j dk ' V ^( k ^ fc '; E ) Xp a (k', k a ) . (2-17) 



Then the momentum representation of Eq. (2-4) reads as 



X 7a (fc, k 



it ) — o 7)Q , 



6(k - k a ) + (-l) Sa V a 5(k + k, 



If we use Eq. (2T8) in Eq. ( |2-17 ) and release the on-shell condition Eq. 
finally obtain 

T 7Q (p, q- E) = V ja (p, q;E)+Y, 7A3 J dk V y/3 (p, k; E) 



(2-18) 



M), we 



/3 

X W kj-l + i e T ^ q ' E) ' 



(2-19) 
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which we call an LS-RGM equation. 

The partial-wave decomposition of the scattering amplitudes etc. are carried 



out in the standard wayll^P . For example, we have 

T ja (k,k';E) = J2 ^T^ sie , <aSe (k,k';E) 
JMM 

xY,^rn'S'S' z \JM)Y eml {k)^rnSS z \JM)Yl 7n {%') , (2-20) 
m' m 

where 7 = [1/2(11) c%, l/2(ll)c2] S' S' Z Y 1 1 ^V' . The prime on the summation symbol 
indicates that £ (£') is limited to such values that satisfy the condition (— 1) = 
(— 1) S V ( (—1)^ = (— l) 5 V ). We write the spin quantum numbers S and S' of the 
T— matrix explicitly for the later convenience, although these are already included in 
the definition of a and 7. The partial-wave decomposition of the LS-RGM equation 
yields 

' ^y]- POO 

T js>e,ase(P, 9; E) = V J s , e se (p, q;E) + ^ — — / k 2 dk 

0f v , K^Y Jo 

x ^s'e',i3S"£"(P^ k'i E) ~^2~ ^2 _ p _|_ i £ Tps»e",ase(k, q; E) , (2-21) 

where E = E™ 1 + (?i 2 /2 / u / g)fc^ and V^ s , e , aSe (p,q; E) is the partial-wave decompo- 
sition of V ia (p,q; E). The Lippmann-Schwinger equation Eq. (2-21) involves a pole 
at k = kp in the Green function. A proper treatment of such siii^ularity for posi- 
tive energies is well known in the field of few-body p roblems .EJ)'|]J)'E2P Here we use 
the technique developed by NoyesB and KowalskiE-§ , and separate the momentum 
region of k (and also p and q) into two pieces. After eliminating the singularity, we 
carry out the integral over < k < kp by the Gauss-Legendre 15-point quadrature 
formula, and the integral over kp < k < 00 by using the Gauss-Legendre 20-point 
quadrature formula through the mapping, k = kp + tan(-7r/4)(l + x). 

2.2. Basic Born kernel 

In this subsection we explain how to obtain the basic Born kernel, Eq. (2-3), in 
detail. In fact the calculation of the kernel is rather tedious and a careful check of 
the calculation must be made. Recent modern techmques p _especially developed in 
the microscopic nuclear cluster theory of light nuclei^' Ef'Eif, have greatly reduced 
the labor of tedious calculations. Since two-body forces used in such applications are 
usually momentum-independent central and spin-orbit forces of the Gaussian radial 
dependence, one needs to extend the technique to more general^wctbody forces 
involving various types of the tensor forces and Yukawa functions EJ)'lH§). Appendix 



*' In this paper we use the phase convention of the S'-matrix defined through the partial wave 
decomposition in terms of the time reversal state i Yi m (r) in the coordinate space. This phase 
convention is different from the standard one by Blatt and BiedenharnliU by the factor , and 
introduces an extra minus sign for the mixing parameter e.j for the TVAL-scattering. Similarly, all 
the Pauli-spinor invariants in Eq. (2-29) are taken to be "real" operatorstJ . 
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A gives a very general formula to calculate the Born kernel directly from such two- 
body forces by simple Gaussian integrations. 

Another important technical development in our quark-model approach to the 
NN and YN interactions is motivated by the rich flavor contents of the spin-flavor 
SUq wave functions of baryons. The operator formalism introduced in Refs.|26|) and 



27) makes it possible to represent this flavor dependence in a transparent form by 
using abstract SU3 operators expressed by the basic electric- and magnetic-type 
SUq unit vectors. This formalism is also useful to deal with the flavor symmetry 
breaking. lJ 

Keeping in mind the operator formalism of spin-flavor factors, we can express 
the basic Born kernel in Eq. (2-3) as 

M B (q f , qi]E) = { e iq f R | G(R, R';E)\ e**'*) 

= M£ N (q f , Qi ) + ME N (q f , m) + Ml N {q f , qi ) S 12 (k, k) 
+ J2M n {q f ,q i )O a (q f ,q i )-e M N (q f , qi ) , (2-22) 

Q 

where e is the relative energy in the final channel (in the prior form) when the channel 
matrix elements are taken at the baryon level. Each component of the Born kernel 
Eq. (2-22) is given in Appendix B in terms of the transferred momentum k = qf — qi 
and the local momentum q = (qj + Qj)/2. In_Eq. (2-22) the space-spin invariants 
qQ = o Q (q f , qi ) are given by O central = 1 andS 

O ls = in ■ S , O ls( ~ } = in ■ S { ~ } , Lsi ~ )a = in ■ S { - ] P a , 
with n = [qi X q f ] , S = -{a x + er 2 ) , = ^(cri - <r 2 ) , 

1 + <Tl • (To , 

and P a = ^ . (2-23) 

For the tensor part, it would be convenient to take three natural operators defined 
by 

O r = S 12 (k, k) , O r = S 12 (q, q) , T " = S 12 (k, q) , (2-24) 

where S\ 2 {a, b) = (3/2)[(<ri-a)(a 2 -b) + (a 2 -a)(ai-b)] — (ai-a 2 )(a-b). The invariant 
Born kernel M fi \qj ,qi) in Eq. (2-22) consists of various types of spin-flavor factors 
and the sp atia l functions f^{0) calculated for the quark-exchange kernel of the 
FB interaction,!***)! It is explicitly given in Appendix B and is generally expressed as 

M n (q f ,q i ) = '£X? f?(9) . (2-25) 
r 

Here we should note that the spin-flavor factors depend on the isospin and that the 
spatial function f^(9) is actually a function of q"j, q\ and the relative angle 9: cos 9 = 

*' The full account of the operator formalism for the spin-flavor factors of the RGM kernel will 
be published elsewhere. 

**) Here we use a slightly different notation for Q from that used in § 2.1, but the correspondence 
is almost apparent. 

***) Here we will show only the expressions for the quark sector. Those for the EMEP sector are 
easily obtained by some trivial modifications. 
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(Jf •Qi r irhe sum over T in Eq. ( [2-25D is with respect to the quark-exchange interaction 
typeset T = E, S, S', D + and Z)_. The factors X|? (possible only for the central 
force) should be replaced with — Xjr,, because of the subtraction of the internal- 
energy contribution in the prior form. Finally the partial-wave decomposition of the 
Born kernel, Vg,p se (qf,qi) can be calculated by usingp 



f: 



Q 

re 



2 J- 



f¥(0)P e (cos9) d(cos( 



(2-26) 



as follows: 

rj cent 
'S'l',S(. 



yj central ft (g/j ft) = £ X ? f? e with ■ <X 2 ) = 2S(S + 1) - 3 



T 



Ys'ffltfafiQi) = $1', its', sSs,i QfQi 



2(2£ + l) 

EyLS f f LS f LS ^ 

A r Ure+i ~ Jre-ii 



irJ LS<-)<t 
v S'l',Si 



Vs'^st^fafiQi) - fis',S&S,l (Si2)i',i 



r 



£ 

r 



+ l) + 2- J(J + 1)] 



x 



( _ 1} 5 2J+1 

(/rj+i " /rj-i) ^ 5, 5 = 1, or 0, 1 , 



Xf s(_) 



LS(-)cr 

r 



„2 T/ T(//) _2ttT(«) 



2£+l 



^+1 + 



(/») 



for 



£±2 and J = £±1 
£ = £' = J J±i 



(2-27) 



Here the tensor component is the sum over fi = T, T' and T" and 



Vi 



T{ff) 



V, 



T(ii) 



x s fse + 



V, 



T(fi) 



x s fse + x s> f, 



1 

9 



T 

s'e 



1 



X S' fs'i + /d+£ + 1 X D- /. 



J i>-e + t; x d- fr>-e 



2X d+ fi 



(2-28) 



Furthermore, (£12)// i is the standard tensor factor. In Eqs. ( 2-25| ), (2-27) and (2-28), 
the spin-flavor factors X^ should be replaced with X^ =1 j- for the EMEP exchange 
terms, fi = CN, SN and TN. Also, the direct term with x = is possible only for 
T = D + . 

2.3. Invariant amplitudes 

Since the three vectors, k = qj — qi, q = (qf + qi)/2 and n = qi x qj = q x k 
are not mutually orthogonal for the YN interaction, we use k = qj — qi, n = q x k 

*' We use the Gauss-Legendre 20-point quadrature formula to carry out the numerical integra- 
tion of Eq. (]2-2e|) over cos0 = -1 ~ 1. 
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and P = k x n = k 2 q — (fc • q)k to define the invariant amplitudes of Eq. ( |2T3| ): 

M(q f , q t \ E) = g + h i(ai + er 2 ) • n + h- i(a\ - <r 2 ) • n 

+h n (<ri • n) ■ (<t 2 -n) + h k (a 1 • k) ■ (<t 2 ■ k) + h P (cr 1 • P) • (<r 2 ■ P) 

+/+ {(o-i • fe)(<7 2 • P) + (o-i • P)(<7 2 • fc)} 

+/_{( < T 1 -fc)(ff 2 -P)-( < Ti.P)(ff 2 -fe)} . (2-29) 

The eight invariant amplitudes, <?0j " ' " > /-; are complex functions of the total energy 
£" and the scattering angle, cos# = g« • qj. Three of the eight invariant amplitudes, 
h-, /+ and /_, do not appear for the NN scattering due to the identity of two 
particles (h- and /_) and the time-reversal invariance (/_ and /+). These three 
terms correspond to the non-central forces characteristic in the YN scatterin g; i.e.. 
h- corresponds to LS^', f + to S^ir^p), and /_ to LS^a, respectively. EZP In 
particular, the antisymmetric LS interactions, LS^ and LS^cr, involve the spin 
change between and 1, together with the transition of the flavor-exchange symmetry 
V 7^ V . In the NN scattering this process is not allowed, since the flavor-exchange 
symmetry is uniquely specified by the conserved isospin: V = (-1) 1_/ . On the 
contrary, these interactions are in general all possible in the YN scattering, which 
gives an intriguing interplay of non-central forces. 

The invariant amplitudes are expressed by the 5-matrix elements through the 
partial- wave decomposition of the scattering amplitude Eq. ( 2-13j ). If we write it as 



M la {q f ,qi;E) = ^ ^ RjS'i',a 



St 



J Mil' 



x^'^'S'S' z \JM)Y e/m/ (q f ) J2(^SS z \JM)Y i * m (q i ) , (2-30) 

rri m 

the partial-wave component R^gij i aSl = (l/2i)(S^ s , e aSe — S^^s' ,s$i' ,i) is obtained 
from the relationship in Eq. ( 2T6| ): 



RyS'i',aSi — - ^ 2-K~h 2 ~ T: r s ' i '> aSe ( q f ,qi,E ^ ( 2 ' 31 ) 

where e 7 = {h 2 /2^)q 2 fl e a = (h 2 /2[i a )q 2 and E = E™ 1 + e 1 = E^ + e a . The 
formulae given in Appendix D are then used to reconstruct the invariant amplitudes 
by the solution of the LS-RGM equation Eq. (2-21). All the scattering observables 
are expressed in terms of these invariant amplitudes. For example, the differential 
cross section and the polarization of the scattered particle are given by 

^- = cr o (0) = | 5o | 2 + \ho + /i-| 2 + \h - M 2 , 
ail 

+ \K\ 2 + \h k \ 2 + \h P \ 2 + |/+ + /-| 2 + |/+ - /-| 2 

P{6) = 2Qm [g (h + h-)* + h n (h - /»_)*] 

+2Sm[h h (f+ -/-)*- h P (U + /_)*] . (2-32) 
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[ k F > £qi case ] 
W(q u q) 



[ k F < £qi case ] 
W(q 1} q) 




k F - £gi k F + £gi 



£qi - k F (qi + k F 



Fig. 1. The weight function W(gi, ?) in Eq. (p-37|) as a function of q. 



2.4. t e '* p prescription for single-particle potentials 

In this subsection we will derive an approximate formula to calculate s.p. po- 
tentials appearing in the lowest-order Brueckner theory. The Bethe-Goldstone (BG) 
equation for the G-matrix solution is given by i3 



' 4:7T f 

G^s'e'^seiPiT'K^) =V 1 J s , e , )aSe (p,q;E) + 7T13 / 

nana ) JO 



pS"i" 



k 2 dk 



j 



Qp{k,K) 



■yS>e',i3S"e" (Pi k] E) J^'^ K-is) ^/ 3S " i "^ aSi ^' ^' ^' u ^ 



(2-33) 



where we assume E = E™ 1 + (h 2 /2p a )q 2 . In Eq. (2-33), Qp(k,K) stands for t 
angle-averaged Pauli operator and eg{k,K;u) is the energy denominator given byo 



efs(k,K;u) = u - E b (kt) - E N (k 2 ) , 
with the s.p. energy E b (k): 

h 2 

E b (k)=M b + —k 2 + U b (k) . 
Zlvl b 

The s.p. potential U b {k) is calculated from 



(2-34) 
(2-35) 



1 PQinax 



jes 



q 2 dq W(q u q) G J aS ^ aS i{q,q;K, 



(2-36) 



where g maa . = (k F + £qi) / (1 + with £ = (Mjv/M a ) and W(<Zi,<z) is the phase space 
factor given by 



W( gi ,q) = ±(l-[-l\x \l) 



X 



W + (i + flV-*£ 

2e(i+o?i« 



(2-37) 



*' For s.p. potentials, we use the notation bi = b and a± — a to specify baryons. 
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with [o|6|c] = max(a, min(6, c))il). The profile of W(qi,q) is illustrated in Fig. 1. 
The starting energy lo in Eq. ( |2-34 ) is a sum of the s.p. energies of two interacting 
baryons: 



lo = E a (qi) + E N {q 2 ) 
= M a + M N + 



K 1 



2(M a + M N ) 



K 1 + 



q 2 + U a {qi) + U N (q 2 



(2-38) 



where K and q are the total and relative momenta corresponding to the initial s.p. 
momenta q\ and q 2 . Once q\ and q are given, the values of K and lo in Eq. (2-36) 
are calculated through 



K -(1 + 6 ql+q 2 -qiq (l + [-l\x \l]) 



(12 



lo = E a (qi) + E N (q 2 ) 



ft? 



(2-39) 



Let us specialize the system to NN and consider the asymptotic behavior of the 
s.p. potential Un(qi) m symmetric nuclear matter. When q\ is sufficiently large, the 
weight function W(q±,q) shows a delta-function like behavior around q = qi/2 with 
respect to the variable q. (See Fig. 1.) We can replace the G-matrix in Eq. (2-36) 
with the one evaluated at q = qi/2, and carry out the integration over q by 



Qmax 



q 2 dq W(qi 



u3 
F 



3 (l + £) ; 



(2-40) 



On the other hand, the G-matrix equation Eq. (2-33) for q — > cxd should approach 
the LS-RGM equation Eq. (2-21), since the Pauli operator Qp{k,K) plays a minor 
role in such a high momentum region, and the s.p. potentials in the energy de- 
nominator ep(k, K;uS) is relatively unimportant in comparison with the large kinetic 
energies. Thus we can expect that the s.p. potential in the high momentum re- 
gion is well approximated by a product of the on-shell T-matrix with q = q\/2 and 
the density of nuclear matter, which is related to the Fermi momentum kp with 
p = (2/3-7T 2 ) kp. It should be noted that only the spin-isospin independent part of 
the T-matrix contributes, since we take spin-isospin sum for the target nucleons. If 
we note the relationship Eq. ( 2-31[ ) and use the partial- wave decomposition of the 
invariant amplitude go given in Eq. (D-2), we can easily show 

U N (q 1 ) = ~U$ r (q 1 ) + ~U} r (q 1 ) with £/*(<&) = 2p g^q = q x /2, 9 = 0) , (2-41) 

except for the overall factor —(2nfi 2 /pq). This relationship is also shown in terms of 
np and pp invariant amplitudes: 

U N (qi) = 0/2) (58(9i/2,0) + So (<7i/2,0)) + (p/2)2^(gi/2,0) 

= (p/2)(g^( qi /2, 0) + g p p ( qi /2, 0)) , (2-42) 
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Table I. Comparison of the phase-shift values and mixing angles (in degree), calculated by the LS- 
RGM formalism (lseq) and the improved variational method (var), for the AN 3 Si- 3 Di - EN 
Si- D\ (upper) and AN 1 Pi- 3 Pi - EN 1 P 1 - 3 P 1 (lower) coupled-channel problems. The mixing 
angles, ei and pi, are only approximately defined for the energies above the EN threshold <-iThe 
phase convention of ei is the standard one for the bar phase shifts. The model is RGM-Ftf. 



Pa 




AN 




EN 




AN 






3 Si 


3 D! 


3 S! 


3 Di 


ei 




lseq 


37.494 


8.656 






13.802 


430 


var 


37.559 


8.654 






13.814 




diff 


-0.065 


0.002 






-0.012 




lseq 


49.169 


15.227 






22.999 


440 


var 


49.263 


15.212 






23.016 




diff 


-0.094 


0.017 






-0.017 




lseq 


12.831 


-5.535 


122.619 


-0.014 


-21.090 


450 


var 


12.825 


-5.528 


122.708 


-0.011 


-21.099 




diff 


0.006 


-0.007 


-0.089 


-0.003 


0.009 




lseq 


14.793 


-4.461 


87.265 


-0.165 


-13.440 


460 


var 


14.801 


-4.458 


87.298 


-0.152 


-13.449 




diff 


-0.008 


-0.003 


-0.033 


-0.013 


0.009 



PA 




AN 






EN 


AN 






'Pi 


3 Pi 


"Pi 


3 Pi 


Pi 




lseq 


9.161 


5.195 






-8.213 


300 


var 


9.161 


5.188 






-8.208 




diff 


0.000 


0.007 






-0.005 




lseq 


82.874 


0.757 






-41.647 


400 


var 


82.605 


0.971 






-41.636 




diff 


0.269 


-0.214 






-0.011 




lseq 


-15.697 


-20.945 


-1.812 


-8.113 


16.559 


500 


var 


-15.690 


-20.944 


-1.812 


-8.120 


16.563 




diff 


-0.007 


-0.001 


0.000 


0.007 


-0.004 




lseq 


-14.077 


-21.357 


-5.725 


-22.826 


13.399 


600 


var 


-14.070 


-21.358 


-5.724 


-22.831 


13.401 




diff 


-0.007 


0.001 


-0.001 


0.005 


-0.002 



where p/2 is the proton and neutron density. This implies that the asymptotic 
behavior of the s.p. potential in the high momentum region is determined by the 
values of the spin-independent central invariant amplitudes at the forward angle. 
This rule is called t e ^ p prescription, which is sometimes used in the relativistic 
mean field theory. For example, if we want to calculate Un(qi) for qi = 6.210fm _1 
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in nuclear matter of the normal density with kp = 1.35 fm 1 (p = (2/3tt 2 ) kp = 
0.1662 fm~ 3 ), we only need to derive the invariant amplitude (<?q(0 = 0) + 3<7q(# = 
0)) for T iab = 800 MeV (q = q t /2 = 3.105 fm" 1 ) and multiply it by the factor, 
-(4-K/q)(h 2 /M N )(p/2) = -(4tt/3.105) • 41.47 • 0.0831 = -13.95. 

§3. Comparison with the improved variational method 

As a check of the LS-RGM formalism, we consider the phase-shift parameters 
of the AN-UN{I = 1/2) coupled-channel system, and compare them with the pre- 
dictions by the improved variational method. Here we consider two different types 
of the AN-UN{I = 1/2) couplings. The first one is the 3 S\- 3 D\ channel coupling 
by the tensor force, which dominantly comes from the EMEP, especially, from the 
one-pion exchange tensor force. The other is the l P\ - 3 P\ coupling by the LS^ 
force originating from the FB interaction. Since the quark model usually predicts 
very strong LS^ force, the Ap scattering observables involve very rich information 
on the characters of the non-central forces in the YN interaction. 

Table I shows a comparison between the phase-shift values predicted by the LS- 
RGM formalism (lseq) and the improved variational method (var) with respect to 
the AN ^S^Di-UN 3 S 1 - 3 D 1 (upper) and AN ^^P^UN l P x - 3 Pi (lower) coupled- 
channel problems. The mixing angle e\ (p\) between the 3 5jand 3 D\ ^Pi and 3 Pi) 
channels of AN is also compared. The model is RGM-FcP, which gives the UN 
threshold energy at pa = 445 MeV/c. In spite of the prominent resonance behavior 
in this energy region, the two methods give very similar values for the phase shifts. 
The difference (diff) is less than 0.1° for the 3 Si- 3 Di coupling. For the 1 Pi 3 P\ 
coupling, we find that the accuracy deteriorates because of the strong resonance in 
the AN 1 P 1 channel. If we avoid this energy region, the accuracy is very good. The 
difference is usually less than 0.01°. 

Since the resonance behavior in the above two coupled-channel systems is largely 
model-dependent, we summarize it in Table II for the three versions of our quark 
model. Here Vzn(i/2)( 3 S) indicates the strength of the central attraction in the 
UN (I = 1/2) channel, which is evaluated from the 3 S effective potential obtained 
through the p = Wigner transform of the exchange kernel. We find that the AN 
3 S± resonance in RGM-F appears as a cusp, when the attraction of the UN (I = 1/2) 
channel is not strongly attractive as in FSS and RGM-H. Similarly, the UN (I = 1/2) 
3 Pi resonance does not move to the ^l./V 1 P\ state in RGM-H, which has the weakest 
central attraction in the UN (I = 1/2) channel among our three models. 



§4. 7V./V invariant amplitudes at intermediate energies 

Figures 2 and 3 compare with the experimental datail) the model predictions 
of the elastic differential cross sections and the polarizations for the the np and pp 
scattering in the Ti a ^ = 400 ~ 800 MeV range. These observables in the lower 
energies are given in Ref.0). The model in Fig. 2 is FSS for the np scattering, and 
that in Fig. 3 is RGM-H for the pp scattering. The Coulomb force for the pp scattering 
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Table II n Resonance behaidor near the EN threshold for the / = 1/2 states, predicted by RGM-FB, 
FSSff'ff and RGM-HEJ. "step" denotes the step-like resonance and "disp" the dispersion-like 
resonance. 



model 


RGM-F 


FSS 


RGM-H 


Ef N (MeV) 

VsN(l/2)( 3 S) 


39 
-38 


77 

-24 


77 

-18 


AN 3 Si 
AN 3 Di 


step 
disp 


cusp 
disp 


cusp 
disp 


EN 3 Si 


S = 180° I 


S < 60° 


5 < 45° 


EN 3 D 1 


s < o 


s < o 


s < o 


AN x Pi 
AN 3 Pi 


step 
disp 


step 
disp 


disp 
disp 


EN x Pi 


s < o 


s < o 


5 ~ 60° 


EN 3 Pi 


8 < 


s < o 


5 < 40° 



is neglected, since it only affects the extreme forward and backward angles in this 
high energy region. The solid curve indicates the predictions obtained by solving the 
LS-RGM equation, while the dashed curve the results of the Born approximation. 
The latter approximation is apparently inappropriate even at these high energies. 
However, it is outstanding that the calculated Born invariant amplitudes, leading 
to these cross sections in the Born approximation, have almost the same order of 
magnitude as the empirical amplitudes determined from the phase-shift analysis 
Note that the polarization vanishes in the Born approximation. In these calculations 
no imaginary potential is introduced. The theory overestimates the np differential 
cross sections at backward angles. The np polarization has an unpleasant oscillation 
around c .m. ~ 110°. There appears a symmetry in the pp scattering because of 
the identity of two protons: The differential cross section becomes symmetric with 
respect to 6 c , m = 90°, and the polarization is symmetric with an opposite sign. 
The pp polarization for Ti a b > 400 MeV shows an oscillation around 90°, which is 
not present in the experiment. Except for these disagreements, the characteristic 
behavior of the energy dependence and the angular distribution is reasonably well 
reproduced within the wide energy range up to 800 MeV. These results indicate that 
the LS-RGM technique is very useful for investigating the baryon-baryon interaction 
above 300 MeV. 

Figure 4 shows the five invariant amplitudes, go(6) ~ hp(6), as a function of 
the cm. angle 6, predicted by FSS for the np scattering at T; a f, = 800 MeV. The 
left column displays the real part, while the right column the imaginary part. The 
predictions by the Paris potential E^P are also shown in the dashed curve for com- 
parison. In these calculations the partial waves up to J = 8 are included. The 



*' Actually Jmax = 8 is not big enough for Ti a b = 800 MeV, as seen from the small ripples of 
the solid and dashed curves in Fig. 4. The partial- wave contributions for J > Jmax from the Born 
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30 60 90 120 150 30 60 90 120 150 180 

6cm. (deg) c .m. (deg) 



Fig. 2. The differential cross sections and polarization for the elastic np scattering at Tj a (, = 400 ~ 
800 MeV. The model is FSS. The solid curve denotes the full calculation with LS-RGM, while 
the dashed curve the one with the Born approximation. Experimental data are from E9) . 



amplitudes are added to obtain the results in Figs. 2 and 3. 
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dotted curve (Arndt) is ±he empirical value, which is calculated from the solution of 
the phase-shift analysis &y by using the real part of the phase-shift parameters up 
to J < 6 and the partial-wave expansion of the invariant amplitudes Eq. (D-2). It 




30 60 90 120 150 

8 c .m. (deg) 



30 60 90 120 150 180 

6cm. (deg) 



Fig. 3. The same as Fig. 2, but for the elastic pp scattering in model RGM-H. The Coulomb force 
is neglected. 
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is clear that the most prominent disagreement between the FSS prediction and the 
other predictions appears in the real part of the spin-independent central invariant 
amplitude Ke go(0) at the forward angle. FSS predicts 3fte g(0) ~ 1.25, while the 
Paris potential and the phase-shift analysis predict Ke <?(0) ~ —1. If we use the t e ^ p 
prescription discussed in § 2.4, these values correspond to the nucleon s.p. potential 
in normal nuclear matter, —17 MeV for FSS and +14 MeV for the latter two, as a 
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Fig. 4. The real (left) and imaginary (right) parts of the five invariant amplitudes, go(9) ~ hp{6), 
for the np elastic scattering at Tj (, = 800 MeV, predicted by model FSS. The dashed curve 
stands for predictions by the Paris potential and the dotted curve those by the phase-shift 
analysis SP82 by Arndt et al. EiP. The partial waves included are J < 6 for SP82, and J < 8 
otherwise. 
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Fig. 4. -continued 
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contribution from the unlike nucleons.til We here find that the attractive behavior of 
the FSS s.p. potentials around this energy region is related to the wrong sign of the 
real part of the spin-independent central invariant amplitude go(6) at the forward 
direction. This difference does not impair the differential cross sections very much, 
since the invariant amplitude has an appreciable magnitude for the imaginary part as 
seen in Fig. 4. This should, however, affect some particular polarization observables, 
and cause a large disagreement compared to experiment. 

§5. Nucleon and hyperon single-particle potentials in nuclear matter 

Figure 5 illustrates s.p. potentials Ub(qi), as a function of the incident momen- 
tum qi, which are obtained for B = N, A and E from the G-matrix calculations 
with the model FSS. The normal density with kp = 1.35 fm _1 is assumed for nu- 
clear matter. The left panels are the results with QTQ prescription for intermediate 
energy spectra, while the right panels are those with the continuous prescription. 113 
The upper panels are for the momentum q\ < 4 fin , while the lower panels for 
qi < 20 fm _1 . In the QTQ prescription, the self-consistency of the s.p. potentials 
is not respected in the momentum region q\ > kp = 1.35 fm _1 . On the other hand, 



*' These are consistent with the numbers given in Table VI. Namely, if we add up the I — 
contribution and one-third of the I — 1 contribution to the real part of the s.p. potentials at 800 
MeV, the potential depth becomes -17.4 MeV for FSS, 15.5 MeV for the Paris potential, and 14.5 
MeV for SP99. 
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Fig. 5. The nucleon and hyperon (A, E) s.p. potentials predicted by the G-matrix calculation of 
model FSS. The results in the left panels are obtained by using QTQ prescription, while those in 
the right panels by the continuous prescription for intermediate energy spectra. The momentum 
interval is < qi < 4 fm _1 in the upper panels, and < qi < 20 fm _1 in the lower panels. The 
Fermi momentum &f = 1-35 fm _1 is assumed. 



it is fully accounted for in the continuous prescription.fcJ This implies that we first 
calculate Un(qi), and then Ua{qi) and Us(qi) are determined self-consistently by 
using the result of Un(qi)- The partial waves included are for J < 9. We can see 
from Fig. 5 that the s.p. potentials predicted by FSS is fairly strongly attractive in 
the momentum interval q\ = 5 ~ 20 fm _1 for all the baryons. In particular, £/jv(<Zi) 
in the continuous prescription becomes almost —80 MeV at q\ = 10 fm~~ . This 
momentum interval corresponds to the incident energy = 500 MeV ~ 8 GeV 
in the NN scattering (see Table III). We therefore need to examine the invariant 
amplitudes carefully in this energy region. 

Let us first examine whether the t e ^ p prescription discussed in § 2.4 is a good 
approximation to the s.p. potentials predicted by the G-matrix calculation. Table 
III shows such a comparison with respect to the nucleon s.p. potential L0v(?i) ( m 
MeV) predicted by the model FSS. Here we find that the maximum value of the total 
angular momentum J max = 7 is actually too small, when the incident momentum 

*•* In Ref. p^[) we have assumed C/s((jri) = (7b(<?i = 3.8 fm _1 ) for q\ > 3.8 fm _1 , in order to avoid 
the unrealistic behavior of the s.p. potentials in the high momentum region. Here qi = 3.8 fm _1 
corresponds to Ti a t = 300 MeV in the NN scattering. 
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Table III. Comparison of s.p. potential Um(ii) (in MeV), obtained by t e ^ p prescsiption (T- 
matrix) and by G-matrix calculation in the continuous choice Ej. The model is FSSIj' : I3'. 

(Real part) 



Tlab 


9i 


T-matrix 


G-matrix (cont.) 


(MeV) 


(fm- 1 ) 


J < 9 


J < 9 


J < 7 


187 


3 


-37.49 


-38.98 


-36.76 


332 


4 


-27.81 


-26.08 


-23.75 


518 


5 


-30.19 


-35.79 


-33.01 


2074 


10 


-79.90 


-78.00 


-54.20 


8295 


20 


-14.81 


-16.07 


-7.82 



(Imaginary part) 

Tiab qi T-matrix G-matrix (cont.) 



(MeV) (fm" 1 ) J < 9 J < 9 J < 7 



187 


3 


-44.91 


-29.71 


-29.53 


332 


4 


-52.88 


-39.64 


-39.43 


518 


5 


-60.51 


-44.68 


-44.88 


2074 


10 


-34.91 


-36.81 


-29.77 


8295 


20 


-22.96 


-21.07 


-2.92 



Table IV. Decomposition of Un(qi) (in MeV) in Table III to the isospin I —M and I — 1 compo- 
nents. The total angular momentum included is J < 9. The model is FSStJ'Q . 

(Real part) 



Tlab 


<?i 


I 


= 


I = 


1 


(MeV) 


(far 1 ) 


T-matrix 


G-matrix 


T-matrix 


G-matrix 


187 


3 


-11.91 


-15.25 


-25.58 


-23.73 


332 


4 


-5.98 


-7.53 


-21.83 


-18.55 


518 


5 


-5.04 


-9.65 


-25.16 


-26.14 


2074 


10 


-28.72 


-27.93 


-51.18 


-50.07 


8295 


20 


-2.37 


-3.60 


-12.44 


-12.47 



(Imaginary part) 



Tlab 


<7i 


/ = 





I = 


1 


(MeV) 


(fm" 1 ) 


T-matrix 


G-matrix 


T-matrix 


G-matrix 


187 


3 


-22.02 


-15.85 


-22.89 


-13.86 


332 


4 


-22.66 


-17.69 


-30.22 


-21.57 


518 


5 


-22.21 


-19.69 


-38.30 


-24.99 


2074 


10 


-15.27 


-14.83 


-19.63 


-21.99 


8295 


20 


-20.62 


-18.78 


-2.34 


-2.29 
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Table V. Reduction factors of the S-meson central attraction due to the momentum-dependent q 2 
term. 





Tlab 


0.0243 q 2 


2 


332 MeV 


0.097 


3.1 


800 MeV 


0.23 


5 


2 GeV 


0.61 


7.5 


4.7 GeV 


1.37 



qi > 10 fm _ . However, if we take the same J max in the two calculations, the accu- 
racy of the t e f f p prescription seems to be fairly good even for the large momentum 
around q\ ~ 10 fm _1 . Quite surprisingly, this approximation is very good even at 
such a small energy as = 200 MeV, as long as the real part of the s.p. potential 
is concerned. This agreement between the two prescriptions becomes clearer, if we 
examine each isospin component with 1 = and 1=1, separately, as seen from 
Table IV. As to the imaginary part of the s.p. potential, the t e ^ p prescription 
seems to overestimate the values by the G-matrix calculation. 

From this comparison and the behavior of the invariant amplitudes in the preced- 
ing section, we have found that the attractive behavior of Un{qi) in the momentum 
interval qi = 5 ~ 20 fm _1 in FSS is related to the wrong sign (of the real part) 
of the spin-independent central invariant amplitude at the forward direction. We 
expect that this situation is common even for the A and E hyperons, and it is a 
flaw of our present quark model (not only FSS, but also RGM-H and RGM-F). 
It should be mentioned that the quantitative aspect of the present non-relativistic 
single-channel calculation in the above momentum region is not entirely trustable. 
The corresponding energy region of the NN scattering is already the relativistic 
energy region, where many inelastic channels are open. Nevertheless, a repulsive 
behavior of the s.p. potential at q\ = 5 ~ 20 fm _1 seems to be a mandatory re- 
quirement even in the single-channel calculation, since channel-coupling effects are 
expected to work attractive to the s.p. potential. Although this energy region may 
already be out of the applicability of our non-relativistic quark model, we need s.p. 
potentials for the intermediate states with quite high momenta when the G-matrix 
calculation is carried out in the continuous prescription. 

In order to solve this problem, we use an advantage of our quark model that the 
effect of the short-range correlation is rather moderate compared with that of the 
standard meson-exchange potentials like the Paris potential. Namely, an improve- 
ment of the Born amplitudes is clearly reflected to the improvement of the solution 
of the LS-RGM equation. In particular, the intermediate-range attraction from the 
scalar-meson exchange has the Born kernel 

q 2 k 2 1 

2M 2+ Md 2 \ (54) 

with k = qf — qi and q = (qf/ + 9i)/2, in the approximation up to the order of (v/c) 2 . 
Here m and M are the meson mass and the baryon mass, respectively. So far we 
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have used only the leading term 1 in the square bracket of Eq. (5-1) and neglected 
q 2 and k 2 terms. A dominant contribution to the spin- independent central invariant 
amplitude at the forward angle comes from the q 2 term, which becomes important 



in the high-energy region. In fact, if we set k 2 
bracket of Eq. (5-1) as 



-m 



as usual, modify the square 



1 



+ 



k 2 



2M 2 8M 2 

,2 



1 



m 



2M 2 8M 2 



m 
8M 2 



1 



1 -m 2 /8M 2 2M 2 



(5-2) 



and redefine g 2 including the (1 — m?/8M 2 ) factor, we find the momentum depen- 
dence like 1 — 0.0243 q 2 for the mass mc 2 = 800 MeV of the e meson in FSS. Here 
\q\ is in units of fm _1 . This non-static term of the S mesons plays a role to reduce 
the strength of the intermediate-range attraction by about 20 % at 800 MeV. (See 
Table V.) Actually the momentum q is not directly related to the total energy, nor 
the direct Born term is good enough to discuss the reduction of the central attrac- 
tion at higher energies. We also have contributions from the inherent zero-point 
oscillation of the cluster wave functions and those from the quark-exchange kernel, 
since our EMEP are acting between quarks. Nevertheless, the discussion here is still 
valid, since the dominant contribution to the intermediate attraction is the direct 
term of the e- meson exchange potential. Bryan-ScottEHP carefully examined these 
q 2 momentum-dependent terms in the S-meson and V-meson exchange potentials. 
They found that the inclusion of these terms has a favorable effect of making the non- 
relativistic approximation uniform to order q 2 , although the main effect is almost 
compensated for by a slight change in the coupling constants. Since these terms ^ixe 
included in the Paris potential^ and all the Nijmegen soft-core potentialsE^'t^P, 
it would be useful to incorporate these terms in our quark model, in order to de- 
scribe correctly the asymptotic behavior of the s.p. potentials in the high-momentum 
region. 

As an example of the quark model with the momentum-dependent q 2 term, we 
show in Table VI the result of the s.p. potential £/jv(<Zi) in the t e ff p prescription, 
calculated by a new model (present). In this new model V-mesons are also incor- 
porated as the EMEP acting between quarks. The phase-shift parameters of the np 
scattering is largely improved in comparison with those of FSS. The details of this 
newmodel will be published elsewhere. In Table VI predictions of the Paris poten- 
tial E2P and of the phase-shift analysis SP99E3, together with the decomposition to 
1 = and 1=1 contributions, are also shown for comparison. We find that the new 
model is still slightly too attractive around T/ a b ~ 800 MeV, but the flaw of FSS with 
too attractive asymptotic behavior in the high-momentum region is clearly removed. 



§6. Summary 



In the quark-model study of the nucleon-nucleon (NN) and hyperon-nucleon 
(YN) interactions, a variational method is usually used to solve an integro-differ- 
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Table VI. Comparison of Un(qi) (in MeV) in the t e ' ' p prescriptkp. and its decomposition ta-the 
isospin 1 = and 1 = 1 components, with other models (ParisEj- 1 ) and experiment (SP99tJ). 
The result by a new model with the momentum-dependent q 2 term (present) is also shown. 



T lab (MeV) 








FSS (J < 8) 






or qi (fm -1 ) 


I = 







I = 


1 


total 






real 




imag 


real 


imag 


real 


imag 


200 


-11.43 




22.08 


-25.47 


-23.48 


-36.90 


-45.56 


400 


-5.11 




22.59 


-23.81 


-33.51 


-28.92 


-56.11 


800 


-5.58 




19.64 


-35.41 


-44.37 


-40.99 


-64.01 


<7i=10 


-17.36 




10.48 


-53.92 


-19.57 


-71.27 


-30.05 


<7i=20 


-7.04 




-1.83 


-12.52 


-2.34 


-19.56 


-4.16 


T lab (MeV) 








present (J < 8) 






or qi (fm -1 ) 


I = 







I = 


1 


total 






real 




imag 


real 


imag 


real 


imag 


200 


-9.40 




16.72 


-26.02 


-20.72 


-35.42 


-37.45 


400 


-1.90 




16.11 


-13.41 


-26.58 


-15.31 


-42.70 


800 


4.18 




15.48 


-4.24 


-38.18 


-0.06 


-53.63 


91=10 


13.44 




11.52 


29.23 


-50.58 


42.67 


-62.10 


gi=20 


39.4 




-47.3 


176 


-150 


215 


-197 


T lab (MeV) 








Paris potential (J < 8) 






or qi (fm -1 ) 


I = 







I = 


1 


total 






real 




imag 


real 


imag 


real 


imag 


200 


-9.23 




17.74 


-27.37 


-19.50 


-36.60 


-37.25 


400 


0.65 




16.87 


-10.76 


-25.09 


-10.11 


-41.96 


800 


10.28 




18.97 


15.68 


-36.98 


25.96 


-55.95 


gi=10 


25.67 




34.83 


78.49 


-89.93 


104 


-125 


gi=20 


32.78 




83.78 


143 


-266 


176 


-350 


T lab (MeV) 






SP99 (Phase shift analysis) (J < 7) 






I = 







I = 


1 


total 






real 




imag 


real 


imag 


real 


imag 


200 


-8.56 




16.37 


-25.68 


-19.76 


-34.24 


-36.13 


400 


3.12 




15.08 


-15.87 


-30.13 


-12.76 


-45.21 


800 


13.28 




16.51 


3.53 


-74.20 


16.82 


-90.71 



ential equation formulated in the (3q)-(3q) resonating-group method (RGM). In this 
paper we have developed a Lippmann-Schwinger formalism forthe (3q)-(3q) RGM 
as an alternative method to the improved variational method c3 . The basic equa- 
tion, which we call the LS-RGM equation, is completely equivalent to the standard 
RGM equation in the coordinate representation on the energy-shell. An advantage 
of solving the RGM equation in the momentum representation is to avoid rapid os- 
cillation of the relative wave functions at higher energies, which leads to the strong 
cancellation of Gaussian trial functions in the improved variational method. This 
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feature of the LS-RGM formalism naturally makes it possible to obtain an accurate 
S-matrix even for the relativistic energies T[ a \, > 350 MeV, where many inelastic 
channels open. Since the S-matrix is very accurate also in the low-energy region, 
an extension to the coupled-channel systems with different threshold energies is very 
successful. 

In this formulation the Born kernel is analytically calculated for all pieces of di- 
rect and exchange terms, which are composed of the kinetic-energy term and various 
pieces of quark-quark (qq) interactions. The qq interactions are further divided into 
the phenomenological confinement potential, a color analogue of the Fermi-Breit 
(FB) interaction, and the effective meson-exchange potentials (EMEP) acting be- 
tween quarks. The Born kernel is then decomposed into partial-wave components 
and the resultant LS-RGMequation is solved by using the techniques developed 
by Noyest!3 and KowalskillS. Since calculations are always carried out in the mo- 
mentum representation, the present formalism has no difficulty to incorporate the 
momentum-dependent qq interaction such as the momentum-dependent Breit retar- 
dation term of the FB interaction, the higher-order terms of the central scalar-meson 
and vector-meson exchange potentials, and the quadratic LS force. A convenient 
transformation formula to derive spatial functions for the direct and exchange kernels 
is given for a very general type of two-body interactions. The numerical evaluation of 
the partial-wave Born kernel is also best suited to convert the LS-RGM equation to 
the G-matrix (Bethe-Goldstone) equation©, in which the Pauli principle is treated 
exactly at the baryon level. 

The accuracy of the LS-RGM formalism has been examined in the AN -EN {I = 
1/2) coupled-channel system. In this system the 3 S±- 3 Di coupling caused by the very 
strong one-pion tensor force yields a prominent cusp structure for the AN phase-shift 
parameters at the UN threshold. On the other hand, a resonance appears either in 
AN Pi state or in EN 3 Pi state by the strong effect of the antisymmetric LS force 
(LS^* 1 force) originating from the FB interaction. The behavior of these resonances 
is rather sensitive to the characteristics of the model, particularly to the strength of 
the central attraction of the EN (I =1/2) channel. We have examined the model 
RGM-F and obtained a satisfactory agreement of the phase-shift parameters between 
the LS-RGM method and the improved variational method, in the wide momentum 
region p tab = ~ 1 GeV/c. 

Using the same parameter set determined at Tua, = ~ 250 MeV, we have 
extended our calculation of the differential cross sections and the polarization for 
the AW scattering to the intermediate energies 2} a & = 400 ~ 800 MeV. Although 
agreement with the experimental data becomes gradually worse as the energy be- 
comes higher, the characteristic behavior of the energy dependence and the angular 
distribution of these observables for the elastic np and pp scattering are reasonably re- 
produced within this energy range. In particular, the invariant amplitudes predicted 
by model FSS for the np scattering at = 800 MeV reproduce reasonably well the 
empirical invariant amplitudes determined from the phase shift analysis, except for 
a few typical disagreements. The most prominent disagreement appears in the real 
part of the spin-independent central invariant amplitude go(9) at the forward angles 
9 ~ 0. This amplitude is related to the single-particle (s.p.) potential obtained from 
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the G-matrix calculation through the t e ^ p prescription, when the incident momen- 
tum is high. The wrong sign of Ke<?o(#) at 8 = in our model is correlated to too 
attractive s.p. potentials in the momentum region q\ = 5 ~ 20 fm _1 . 

We have also examined the accuracy of the t e f f p prescription by using the G- 
matrix solution c3 and the present LS-RGM formalism. This prescription is a good 
approximation for the real part of the s.p. potentials in the energy region from 
Ti a b = 200 MeV to several GeV (q\ = 3 fm _1 ~ 20 fin -1 ), as long as the maximum 
value of the angular-momentum cut for partial waves is commonly taken. On the 
other hand, the imaginary parts of the s.p. potentials are usually overestimated in 
the t e f ' p prescription. 

Perhaps the most striking feature of the quark model found in this investigation 
is a moderate effect of the short-range correlation. In the standard meson-exchange 
models, the observed phase-shifts are reproduced as a cancellation of very strong 
repulsive and attractive local potentials, and the Born amplitudes of the Paris po- 
tential, for example, are one or two order of magnitude large, compared to the 
empirical invariant amplitudes. On the other hand, the short-range repulsion in the 
quark model originates mainly from the nonlocal kernel of the color-magnetic term 
of the FB interaction. Born amplitudes of the quark model therefore have almost 
the same order of magnitude as the empirical amplitudes obtained by solving the 
LS-RGM equation. This implies that the short-range correlation in the quark model 
is rather moderate compared with the meson-exchange models. It can also be seen 
in our recent quark-model study of the s.p. spin-orbit potentials for the nucleon and 
hyperons. In Ref.|35|) we have calculated the strength factor, Sb, of the s.p. spin- 
orbit potentials by using the G-matrix solutions and found that Sn does not obtain 
much effect of the short-range correlation, on the contrary to the standard potential 
models like the Reid soft-core potential with the strong short-range repulsive core. 
Since the Born amplitudes in the quark model reflect rather faithfully characteristics 
of the LS-RGM solution, it is easy to find missing ingredients that impair the model. 
In fact, we have discussed that the wrong sign of the invariant amplitude $te go(9) in 
our quark model is related to our neglect of higher-order momentum-dependent cen- 
tral term of the scalar-meson exchange EMEP. A preliminary result of a new model, 
which incorporates this term as well as the vector-meson EMEP, shows that the s.p. 
potentials in the t e f f p prescription have a correct repulsive behavior in the asymp- 
totic momentum region. The details of this model will be given in a forthcoming 
paper. 

Appendix A 

A transformation formula to the Born kernel for momentum- dependent 

two-body interactions 

As to the general procedure how to calculate the RGM kernel for two-cluster 
systems composed of s-shell clusters, Refs. and |2^) should be referred to. Here 
we use the same notation as Ref. 23) and give a convenient formula to calculate the 



Born kernel for two-body interactions with momentum dependence. 

Suppose a two-body interaction of the RGM Hamiltonian is given by Vij = UijWij 



Lippmann-Schwinger RGM Formalism 



27 



with the spatial part u = u(r,d/dr) and the spin-flavor-color part Wy = wfj w£. 
The full Born kernel for this interaction is defined by 

6 

(e^ 4>\J2vi^'\e iq ^ r <t>) , (A-l) 

i<j 

where 4> = (j) s P ace £ with £ = £ £ is the harmonic-oscillator (h.o.) internal wave 
function of the (3q)-(3q) system and A' is the antisymmetrization operator between 



two clusters. This expression is reduced to the form Eq. ( 2-14 ) by the use of the 
double coset expansion A' — > (1/2)(1 — 9P3e)(l — Pq) with Pq = P14.P25-P36) yielding 
the basic Born kernel Eqs. (2-22) and ( |2-25| ) with 

M(q f ,q l )=Y J X xr M xT {q f ,qi) . (A-2) 
xT 

Here the superscript 17 for specifying the type of the interaction is omitted for 
simplicity, and the sum over x is only for x = and 1. Furthermore, X x t is the 
spin-flavor-color factor defined through 

r 

X x t = C x (z x £ | ^2 w ij I 

i<j 



for x= , (A-3) 



where x = with zq = 1 and Co = 1 corresponds to the direct term and x = 1 with 
Z\ = P36 and Ci = —9 the (one-quark) exchange term. The suffix T stands for the 
interaction type T = E, S, S' , D + or D_, for which the pairs are properly 
selected. The color-factors X^ r = (z x £ c \ \ (p) for each G T are given as 
follows. For the quark sector with wfj = (l/A)(XfXf), X$ E = -(2/3), X^ D+ = 0, 
and X? E = Xg = Xg, = -(2/9), X? D+ = 1/9, X? D _ = 4/9. For the EMEP 
sector with wfj = 1, X^j- = 1 and Xyj- = 1/3. The spin-flavor-color factor for the 
exchange normalization kernel is defined by Xn = (— 3)(P36 £, SF \£, SF )■ We also 
need X K = 24(P 36 £ 5F | Y(6) - Y(5) | £ SF ) for the exchange kinetic -energy kernel of 
the YN systems. On the other hand, the spatial part of the Born kernel is defined 
by 

M xT (q f ,qi) = (z x e i «rr (j >s P a C e\ u ..\ e i qi .r (l) s P <Ke ) ^ (* J) G T , (A-4) 

which we now evaluate. 

The standard procedure is to use the h.o. generating function 



A,(r,z)=(fiy e -7(-*/v^+* 2 /2 ; (A . 5) 

and first to calculate the so-called complex GCM kernel defined by 

I xT (z, z') = { z x Ay(r, z) ^ ace I ua I A 7 (r, *') ^ pace ) . (A-6) 
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This GCM kernel is divided into the general form of the norm kernel and the inter- 
action function T(z,z'): 

I xr {z,z') = e^-f> z - z ' T(z,z') , 
T(z,z>) = <(0*)flj(0*)fl£ )- 1 ((0 S )5 7 |(0 S ) 5 p- 1 

x((Os)sJO S )sJu|(Os) 5 ,(0 S ) 5 ,> . (A-7) 

Here (Os)s stands for a localized (Os) wave function around x = S, and S a etc. 
are the generator coordinates specifying the position of clusters c3> . For a simple 
Gaussian two-body interaction u(r) ~ e~ Kr , T(z,z') is given by 

3 

T(z,z') = ( ^ 2 e~^ ( P z * +qz '^ x (polynomial terms) , (A-8) 

where v = 1/26 2 is the h.o. constant of the (Os) clusters and A is given by A = 
(1/2//)(k/ (v + k)) with 7 = fiv and /u = 3 • 3/(3 + 3) = 3/2. The parameters p and q 
are 0, 1 or —1, depending on T, which are explicitly given in TABLE II of Ref . |23|) . 
The transformation to the Born kernel is achieved through 

3 

M xT (Qf,Qi) = 2 e^ {q2 f +q2) J da db e -i(" 2 + &2 ) ^ir^m-b 

xl xT {^a,^b) ■ (A-9) 

When u involves momentum-dependence, it is convenient to use the momentum 
representation and write 

(P\u\p') = p#',9') with k' = p-p', q' = ~(p + p') . (A-10) 

Then the interaction function in Eq. (A-7) is given by 

T (z,z') = JL f dk ' dq> e-Ho' 2 +W 2 ) e \^ + i q '-r e «.x u(fc / ^ 

(27T)' 3 \7TV J J 

with r = — l —{pz* - qz'), X = --^—(pz* + qz') . (A-ll) 
y/l 2^7 

If we use Eq. (A-ll) in Eq. (A-9), the integration over a and b can be carried out 
and we obtain the following formula after some rearrangement of terms: 

**■(«,•) = «."<«,.) (i)7*'^{" (' + 1) i"' 2 " ^™'} 
i i \ 3/2 r . , r i i 



-irv 1 — a/2/j, 
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Here M% (qf,qi) is the normalization kernel given by 

with t = 1 — x / ft, k = qf — qi and q = (qj + qfj)/2, and the coefficients appearing 
in Eq. (A- 12) are defined by 

p 2 + q 2 — 2rpq „ p 2 + f/ 2 + 2rpg p 2 — q 2 

a = » , a = = , £ = 



1-r 2 ' 1-r 2 ' 1-r 2 ' 

Let us specialize the two-body interaction to the Yukawa function and the Gaus- 
sian function: 

u(k,q) 



k 2 + r" 2 

J_ 1,2 



u(k,q) = { K ^ m . (A-15) 



{lYe~- K k u(k,q) 

Here u(k, q) stands for a polynomial function of k and q, and the degree of the 
polynomial is usually at most the second order in q. The Born kernel for the Yukawa 
function is calculated from the formula for the Gaussian function through the integral 
representation 

3 9 

Air 1 f°° . _j_ 1 7T \ 2 _j_ fc 2 ... / mu x 



m^= / du e ^? - e 4« with k = — . (A-16) 



fc 2 + m 2 y 7 ^" 7o 

The resultant Born kernel is expressed by generalized Dawson's integrals h n {x) and 
modified Yukawa functions y a (x), Z a (x) etc., which are given by the error function 
of the imaginary argument: 

h n {x) = (2n + 1) e~ x2 f 1 e xH2 t 2n dt = (2n + 1) ^— F e u ' u 2n du , 

JO x zn+i J 

y a (x) = e a ~* [\-% +xH2 dt , 



t ( x ) = e a ~ x2 F e -$ +x2t2 t A dt . (A-17) 
Jo 

For the Gaussian kernel Eq. (A- 12) is reduced into 

3 

^ (w , tt )_^ (w , gl )(-i- T J 3ig )V{i T ^=v«} *><„,«) , 

(A-18) 

where the polynomial part is given by 
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1 1 V /2 ( J ' f 1 1 fll 
— i dqexp< — q \ 



xu Ik' + 2^yV,q' - — k' + 



4/i 2^7 

Here q') = u(—k', —q 1 ) and the simplified notation 

A 



1 + Xa 



V , W = A-eV 



(A-19) 



(A-20) 



are used. If u(k,q) does not involve g-dependence, the q' integral is carried out in 
Eq. (A-19) and we obtain 



V{q f ,qi) 



1 + Xa\2 
87T7A 



J k' exp|-i±^ fc' 2 } u(k' + 2^v) . (A-21) 



Appendix B 

- Born kernel - 



In this appendix we show the invariant Born kernel Eq. ( |2-25 ) obtained by us- 
ing the transformation formula in Appendix A. These are functions of k 2 , q 2 and 
k • q with k = qj — qi, q = (qj + qi)/2, and expressed by the special functions 
given in Eq. (A-17). Only 6 with cos# = qf ■ q% is explicitly written in the spatial 
functions f-x(O). The S'-type spatial function f§r(0) is obtained from f§{9) with 
k — » —k. There is no E-type possible for the non-central forces. The partial wave 
decomposition of the Born kernel is carried out numerically through Eq. ( 2-26| ). The 
spin-flavor-color factors in the quark sector are also shown in the operator form in 
the isospin space. 



B.l. Quark sector 

Exchange normalization kernel 



M N (q f , qi )=X N f(9) , 



with 



/(fl) = (v&) 3 exp l--{q 2 + k 2 ) 



(B-l) 
(B-2) 



Exchange kinetic-energy kernel 



M K (q f , qi )=X N 



Y 



f Kl (9) + X K -(l 



where A = (m s /m u d), Y is the hypercharge of the total system, and 



-x 2 m ud c 2 f(6) 



-l + lb 2 (2q 2 + k 2 ] 
1 + U 2 (q 2 - k 2 ) 



(B-3) 
(B-4) 
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with x = (h/m U( icb). 
Color-Coulombic term 
M cc (q f , qi ) = X N 

where 



2 [fi c (0) - f s c {e) - f s P{e)) + f D c + {0) + fE c (e) 



(B-5) 



2 asxm ud c 2 ^-(V3irb) 3 

7T 6 



X < 



(&) 5 exp {-^ 2 [f (g 2 + fc 2 ) - fc • g] } Kq (^6|g + fc|) 

(i)%xp{-l6 2 ( g 2 + ife 2 )} ^o(^|) forT 
(l)%xp{-i6 2 fc 2 } /io(^%|) 



The £-type factor is given by /f c (6>) = v/2/Va5xm Md c 2 (4/3) /(0). 
Breit retardation term 

M^(q f , «) = If (0) - |/f C + - 1/. 



5 

D+ 
D- 

(B-6) 



E 



+X MC f MC {9) + X MV f MU {&) + x «0 ^C-^) 



MC ( MC/ 



-MC ;MC/ 



where 



IT 



(6) = ^l as x 3 m ud c 2 ^{V^b) 3 
exp{-±6 2 (q 2 + fc 2 )} [§-(6q) 2 ] for T 
(fT)"exp{-^6 2 [§( g 2 + fc 2 )-fc- g ]} 

i7 (^b\q + fe|) - f b 2 (k ■ q)^ (-^b\q + k\) 



(0) 



(B-7) 



x < 



+ ^b 4 (q 2 + k-q)(k 2 + k-q)H 2 (^ n b\q + k\)} for T = S 



(i)%xp{-l6 2 r g 2 + Ifc 2 )} 
+ l(bq) 2 Hl (±6|fc| 

(|)%xp{-I^ 2 } 



(|6|fe| 
6 4 (fc.q) 2 # 2 (i&|fc| 



(B-8) 



for T = D A 



2H Ub\q\)+^bk) 2 H 1 (^b\q\) 



for T = D. 



\ b \k.q) 2 H 2 (^b\ q \) 

The functions H n (x) with n = 0, 1, 2 are expressed by /i n (x) through 

Hq{x) = ho{x) 
Hi(x) = ho(x) + - hi(x) 

H 2 (x) = ^hi{x) - ^h 2 (x) . 



(B-9) 
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The function /f c (9) is given below. 
Color-Magnetic term 

M GC (q f , qi ) = -X^f^{9) + £ X<fCfGC {9) 

where /f c (6) = ^/2/TTa s x 3 m ud c 2 f(0) and 
fGC, 



S 

D+ 



fr(0) = \ -a s x 6 m ud c 2 (V3TTb) 



x < 



TT ex P 



4/„2 
3 l 



<Z 2 + fc 2 ) - fc • q] } 
(I) f exp{-i6 2 (g 2 + ifc 2 )} forT 
(§) f exp{-±W} 



(B-10) 



(B-ll) 



LS term 

In general there are three different types of LS terms in the YN interaction; 
i.e., ft = LS, LS(~^ and LS^cfcZP These are different only for the spin-flavor-color 
factors. For each type we have 



M n (q f , qi )= ]T X?# s (9) , 



where 



-a s x m ud c -(37r)2 6° 

7T Z 



X < 



(A) 5 exp {-£& 2 [f (g 2 + fc 2 ) - fe - q] } ~h± (^b\q + fc 
(i) f exp{-l6 2 (g 2 + Ifc 2 )} ^(I6|fc| 
(l) § exp{-iW} ^(^1 



(B-12) 



for T 



S 

D+ 
D- 

(B-13) 



Tensor term 

There are three different types of tensor terms as in Eq. (g2g).ED These 
given by 



arc 



x^^(e) + x 1 s ,f 1 S r(e) + i 



M T (q f ,q t ) 1 

M T '(q f , qi ) J-"»^v^ ■ -a^v^ ■ | v / //( (#) 



(B-14) 
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where 

/?(*) = -\[^asx 3 m ud A37T)h 5 

V 7T 5 



X < 



(n) 5 TT expj-^^ 2 [§( Q 2 + fc 2 ) - k ■ q] } h 2 (^b\q + k\) 

(i) § iexp{-l6 2 ( g 2 + ife 2 )} forT 



5 

(B-15) 



B.2. EMEP sector 
Direct terms 

The S-meson central (CN), PS-meson spin-spin (SiV) and PS-meson tensor 
(TN) terms are given by 

Mg N ( qf , qi )=X™fE N (9) , 
M S D N (q f , ft) = X D+ (tT! • <X 2 ) /£» , 



(B-16) 



where 



fg N (9) = -iirhc- 



m- 



f s D N {9) = AttHc^ 
fn N (0) = ~^hc 



+ k 2 
1 



1 

3 \m 2 + k 2 
1 1 



3m 2 rn 2 + fc 2 



fE N (0) 



3m 2 



(B-17) 



Here m (= mc/fi) is the meson mass, and the relationship Xq$ + = Xq D+ (cti • er 2 ) is 
used in Eq. (B-16). The factor c$ in Eq. (B-17) is a reduction factor for the spin-spin 
contact term of the PS-mesons. 



Exchange terms 



M( C s%\q f , qi ) 

M TN (q f ,q t ) \ 
M TN '(q f ,q t ) j 

rTN" 



-\ civ > r\ SN > l 



xi s f 1 s N {e) + x{ s ,f^{9) + 



M TJS " (q f , q t ) = 2 Xf s f™(9) - X( s , (9) 



rT fTN , 



where 



-47T 



3\! 



hcb 2 



(B-18) 
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x < 



exp{-i&V + fc 2 )} 



£)% X p{-£& 2 [§(g 2 + fc 2 )-fc.q 



for E 







i 



(|)%xp{-^(^ + ife 2 )} 5U + (i6|fc|) 
(|)%xp{-l6^} 



C <5 2q Dj 



for T = S 
for T = D+ 



7^1 



2a r 



for T = L>_ 



3\ 2 ^ ,., 



-4tt ( - hcb z — 
2 J \mc 



x < 



3 

JL)- 2 ±exp{-±b 2 [l(q 2 + k 2 )-k-q}} Z as {-^b\q + k\) for T = S 
I) § iexp{-I6 2 ( q 2 + Ifc 2 )} ^ D+ (j6|fc|) for T = D+ (B-19) 
(|) f Iexp{-iW} Z aD _(^6|g|) for T = D_ . 



The coefficients a-r are given by as = as' = (11/8)oe, od + = 2a E , and a E _ = 
(3/2)a E , with a E = (m&) 2 /2 = (l/2)(mc&/ft) 2 . 

B.3. Spin-flavor-color factors in the quark sector 

The spin-flavor-color factors in the quark sector for the NN system are given in 
Eq. (B.6) of Ref.|5|). We should note that X r there is actually X^ c for T = D+ 
and D_, and = X§, c = -X E . (Note that X E is not X^ c .) Furthermore, 



XM C = XM C = X N , X S IC = AfF = -2X N , X% IC = -X", and X K = for 



NN. The spin-orbit and tensor factors for YN systems are found in Appendix C of 
Ref. 17). Here we list up the one-quark exchange (x = 1) central factors {X^)b 3 Bi 
in Eq. (A-3), with respect to the quark sector of the AN, UN and SN systems. The 
spin-flavor factors in the EMEP sector and the detailed derivation of these factors 
will be published elsewhere. In the following, the parameter A (= m s /m U( fl) controls 
the flavor symmetry breaking of the FB interaction. The spin and isospin operators, 
<Tj and Tj, are with respect to the two baryons Bi with i = 1 or 2 in the initial state, 
and the flavor exchange operator Pp is supposed to operate on the ket state. We 
assume B 2 = B4 = N. The spin- flavor-color factor for the exchange normalization 
kernel is given by X^ = X D f_(X = 1). The factors of the interaction type T = S' 
are discussed in S 5 of Ref. IT7T) . 



[ B 3 B 1 = AA } 



X K = l-(l + a 1 -a 2 )P F , 



X 



MC 
D- 



2+-{l+a 1 -a 2 ) P F 



*•* A = m u d/m a in Ref. p"7|) is a misprint. 
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Y MC 



Y GC 



Y GC 



X 



GC 



1 + J + (1 + <Ti • <T 2 ) Pf 
3 + j + U + (1 + O"! • <T 2 ) P F 



1 1 

2 + 24 
1 

72 



4 1 \ , \ ( 4 1 



1 

"72 



( 13 + J2) + + • °"2) + 6Pf) 



^ + 3 ( 3 + ^" 



(B-20) 



[ B 3 B 1 = EE 



11 2 
X K = 1 + -n ■ r 2 - -(1 + ri ■ t 2 )Pf + (o-i • cr 2 ) - 



5 1, 

1 + -Ti ■ T 2 - -(1 + Ti • T 2 )Pf 



Y MC 



2 
9 

4v 1 + i 



+ (<Ti • er 2 ) 
1 



i + ^n •T 2 + T ^(i + ri -T 2 )P F 

, 1 / A 2 A 5 

1 + -Ti ■ T 2 + (O-i ■ <T 2 ) - I 1 + -Ti ■ T 2 



"12 ( 1 + ^ CTl '' T2 ) (1 + Tl ' T2)PF ' 



1 . ,2/5 

1 + gTi • T 2 + (CT! • CT 2 ) - I 1 + -Ti • T 2 



+ T2 (* + 1) ( X + S* 1 ' °" 2 ) (1 + Tl ' T2)PF ' 

^ = 1(1 + ^.^) +1(1 + 1 + ^) (l+n-T^ 



x 1 

+ (O-i • <7 2 ) — 



5 1 / 4 1 > 

l+3T 1 -T 2 + -(l-- + xf )(l+T 1 .T 2 )P F 



D+ 3-72 



23+ A + Fj + 3 131 



1 



8 9 \ 34 /, 

X + A2 ) Tl ' T2 + y (1 + Tl ' T2)Pp 



2 3 \ 1 / 2 15\ 11 



{ ( l7 - f + ^) + C 11 - ^ + ^; ^ • " 2 

1 



2 1 



1 / 8 3 \ , . „ . . 1 

+2 ( 39 -^ + ^) t '^ + K i -M)< i+ti - t2 H} • (B ' 21) 
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[B 3 B! 



X K = 1 + in • r 2 - (<ti • <r 2 ) i ^1 + jjri • r 2 ^ , Xat = -ix x , 



MC 



MC 



- x x N , 



Y MC 

A s 



1 + j ) X N 



vGC 
A D_ 



1 + qTi ■ T 2 - (<Ti ■ CT 2 ) — ( 1 + -Ti ■ T 2 



GC 



1 



3 • 72 



n 8 9 

9+ A + A^ 



+ 



9 -A + A^ 



- (<Tl • <T 2 ) - 



A A 2 



n • n 

8 15 



3-T + T^ +7T 1 5 + - + - 7 n-T 2 



A A 2 



A S 



72 



A A 2 



^ 5-1^ +13- ^r,, 



A A 2 



A A 2 



40 15 



13 — - + —77 ) + — ( 41 — - — |- — r I n • n 



A A 2 



[ = /lZ 1 and 17/1 ] 

X K = P F - (<n • a 2 ) i (J + Pf) , 



A D_ 



A S 



A D_ 



1 

4 

_5_ 

36 



-P F + ((Ti • <T 2/ 



1 /10 1 



3"A P - 



1 + 1) ((Ti • o- 2 ) + i (l - i(Ti • er 2 ) P F 



5 

36 



3 + x )((t 1 -(t 2 ; 



1 + 



A 



GC 
0+ 



2 • 81 
1 

3-36 



•Ti • <X 2 

2 
A 



1 

24 



n 4 1 
1 + 



1 - • <t 2 ) p F , 



- (<n • o- 2 ) - 1 - — + 



(B-22) 



1 



3A A 2 



+ 9P F ) + (o-i • o- 2 ) I (49 + i| - 18P F 



All 



+ 9-4 + 



Pf 



T2)J 



+ (<7i-<7i)J (39- g + 



A 1 



7 



Sid 

A- + A* 



Pf 




Pf , 



(B-23) 
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Appendix C 

Partial-wave decomposition of the Born kernel ■ 



In this appendix we derive some formulae for the partial-wave decomposition of 
the Born kernel. Suppose the basic Born kernel Eq. (2-22) is expanded into partial 

waves as 



M(q f , qi ) = £ M Q (q f , Qi ) Q {q f , qi ) 

a 

= AlT R S'£'Mlf^i) y(i'S')JM{qf,spin)y^ s) j M (qi;spin) , (C-l) 



JMU'SS' 



where yus)JAl(Q'i spin) = [Yi(q)xs(spin)] jm is the standard space-spin function. 
We can show that the partial-wave components are given by 



RJ S'l',sMjiH) = ^2 R S'i',Si(.1f' < li) > 



R s> J e>,se( < lfiq.i 

with cos = (qf ■ qi) and 
4ir 



2 J- 



M S2 (q f ,q i )I s J4, Si (e)d(cos9) 



(C-2) 



tsd'M ) = 9 7,1 J2{y(£'S')JM(Qf, spin)\0 Q {q f , qi)\y^s)JM{Qu spin)) spin , 
M 

(C-3) 

by using the rotational invariance of Rg/ji seilfi Qi)- I n Eq. ( C-3 ) the matrix elements 
are with respect to the spin variables. For the spin-independent central (O = 1) 
and the spin-spin (O ss = (tri • (T2)) terms, the factors Eq. (C-3) are given by 



h>tr&W) = &s>,s6vA [25(5+1) -3] f Pl{ 



cos ( 



(C-4) 



For the LS and tensor terms with Eqs. (2-23) and ( |2-24| ), some calculations using the 
angular-momentum algebra yield 



Is't',Si{0) = -^S,l^S',S^£',£QfQi2 



1 + 



2- J(J + 1) 



+ 1 



r ls(-) \ 

\LS(-)af- 



x sin 9 P/(cos 6) for £ = J, J ± 1 



mi 



e,st ( d ) ~ s J,e s £',£ I [5(5 + 1) - 1] 

for 5', 5 = 0, 1 or 1,0 , 

I \ (Si2)ie [g/ 2 P/(coe0) + qi 2 P e/ ( C os9) 



sin 6 Pj (cos 



f Pj (cos 6) 

T^QfQi j cos 0p e ( cos 0) + |j^jySU10P/(cOB0) 
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for 



£ ± 2 and J = £ ± 1 
£ = £' = J, J ± 1 



g/ Pi(cosO)-qi />(cosi 



(C-5) 



Here Pj (cos 1 



; sin^Pj(cos^) is the associated Legendre function of the first kind 
with degree 1, and (S12)/' « is the standard tensor matrix elements. The final braces 

in I\ei \o (9) in Eq. (C-5) are applied to the choice of £ and I . Note that the LS 
and tensor components are non-zero only for the spin triplet states, while the LS^~' 
and LS^~'(T components are accompanied with the spin transition between S = 
and 1. 

The partial-wave components of the full Born kernel, Vg,g, se (qf, qi) in Eqs. (2-27) 
and (2-28), are easily obtained from the above formulae. 



Appendix D 

Partial-wave decomposition of invariant amplitudes 



Reconstruction of the invariant amplitudes in terms of the 5-matrices is an 
opposite process to the partial- wave decomposition of the Born amplitudes, given 
in the preceding appendix. We can use the same formulae, Eqs. ( |C-4| ) and (C-5), 
to derive the following result for the partial-wave decomposition of the invariant 
amplitudes. 

Suppose the invariant amplitudes Eq. (2-29) in the isospin basis has the following 
decomposition: 



McaiQf, Qi) = V (1 + 5 C1 , C2 )(1 + S ai,a 2 ) 

x 47r R yS'£',asey(e'S')JM(Qf-,spin)y* es)JM (^spin) . (D-l) 

JMU'SS' 

Here a = (01,02) and c = (01,02) imply two-baryon configurations corresponding 



to the initial and final states, a and 7, respectively. The front factor, J (1 + S Cl ,c 2 ) 



x y (1 + 0ai,a 2 )> is a factor for two identical particles, and is 2 for NN, for example. 
Without this trivial factor, the invariant amplitudes for a fixed ca are given by 



9o 

9ss 



jes 13 



1 



1 (2J + 1) 

4^ J(J + 1) 



3 [25(5 + 1) -3] 

RlJAJ P}(COS0) 



R J si,siPi(cc*e) 



/?._ 

f- 



JR J u+i,u+i p ) + i( cos ^)-( J + 1 )^J-i,i^i p }-i(c°s0) 

Pj (cos 9) , 



ly, (2J+1 



vW+Tj 

l^p (2J + 1) 



R(j,oj 



Rqjaj 



Rij,oj 



Rqjaj 



P}(cosi 
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h ( k 



hn — 9ss 2 



cos 9 
sin 9 



T <^ cos(6» - 20) 



sin i 
cos 9 



sin 6* sin 9 

COS9 /»(!) + J- /»(3) 



sin(6> - 20) / 



(i) 



sin# 



/+ = sin(# - 20) 



cost/ 
sin# 



sin# 



1 



sm( 



+ cos(9- 20) /W, (D-2) 



with 



^ (0) = ^E [(2^ + 1) R J u,u Pj(cos9) - (J + 2) R( J+1M+1 Pj + i(cos9) 



-(J-l) R(j- hlJ -i Pj-i(cos9) 



E 



(2J + 1) 



_J(J + 1) 



1 /?J 



1J+] I ./ 



-1 + #1J- 



1,1.7+1 



4^^7(7+1) 

ly, (2J+1) 
/I 2^ 



R 



j 

u+i,u- 



.i + R 



P'(j+i,u-i 



u-i.ij+i 



^ij-i,u+i 



^ij+i,u+i -fj+i( cos( 



Pj(cos0) , 
Pj(cos0) , 

P](cos#) . (D-3) 



In Eqs. (D-2) and (D-3) we have used shorthand notation R J S i e se to denote R^g'v aSi , 
where the spin values S' and S are uniquely specified from 7 and a, respectively. In 
Eq. (D-l) the summation over V' and V is already taken by adding up over all possible 
values of £ and S, for which V is uniquely determined through the generalized Pauli 
principle, (— l) e (— 1) 1 ~ S V = —1. Note that a simple relationship, hj,+hp+h n = 3g ss , 
holds. The on-shell scattering amplitudes are conveniently expressed by only two 
angles 9 and 0, the latter of which is defined through cos = P • qi with 



cot 



sm ( 



\qi\/\qf \ -cos9 



(D-4) 



For the elastic scattering with ai = c\ and ai = C2, the magnitude of the cm. 
momentum is equal before and after the scattering (|<jfj| = |q/|), thus yielding = 
9/2. In this case, the time-reversal invariance gives a further simplification for non- 
zero invariant amplitudes. Since the 5-matrix is symmetric, /_ in Eq. (D-2) and 
in Eq. (D-3) are both zero. Combining these two conditions, we find that /_ = /+ = 
for the elastic scattering. This is the well-known result of six independent invariant 
amplitudes for the elastic scattering with the time-reversal invariance. 
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